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In studying the formation of topological defects, it is conventional to assume the "geodesic rule" which is equivalent to mini- 
mizing gradients of the order parameter. This assumption has been called into question in field-theoretic studies of first order 
phase transitions and in the case of local defects. We present a scheme for numerically investigating the formation of strings 
without assuming the geodesic rule. Our results show that the fraction of string in infinite strings grows as we deviate from 
the geodesic rule. 
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Within all known theoretical frameworks, the early 
universe must have been a close analog of condensed mat- 
ter systems that we now study in the laboratory. In par- 
ticular, the early universe must have undergone phase 
transitions during which it is very conceivable that topo- 
logical defects would have been produced. The observa- 
tion of any left-over defects today would directly provide 
us with a window to the universe when it was a mere 
fraction of a nanosecond old. 

The study of topological defects has made long strides 
over the last twenty years [Q. Starting with order-of- 
magnitude estimates, we are now in a position to com- 
pare theory with laboratory experiments and astronom- 
ical observations. As our tools get sharper, finer details 
of our understanding of the properties of defects will be 
tested. 

The properties of defects at formation, i.e. at the 
phase transition, are an important ingredient in their 
cosmology. For example, if the location of magnetic 
monopoles is very strongly correlated with the location 
of antimonopoles, rapid annihilation would ensue and the 
monopole over-abundance problem in cosmology would 
be absent. Similarly, the presence of infinite strings in 
cosmic string networks is expected to be very important 
in the cosmological role that they can play. 

The properties of defects at formation have mostly 
been studied via lattice-based numerical simulations us- 
ing the so-called Vachaspati- Vilenkin (VV) algorithm . 
To describe these efforts on a concrete footing, consider 
the global U{1) model 
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where (j) = |(/)|e*" is a complex scalar field. Apart from 
the discretization of space in the form of a specific lattice, 
the simulations also assume: 

• Only the dynamics of the phase of <j) is relevant. In 
other words, is set to be rj everywhere. 

• The "Geodesic Rule" : The gradients of a cost a lot 
of energy and are minimized. 



Then the simulation proceeds by laying down phases a at 
every lattice site, finding phase differences by minimizing 
gradients and then evaluating 
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around all plaquettes of the lattice. The value of this inte- 
gral yields the winding number and hence the number of 
strings passing through each plaquette. Then the strings 
are joined together and in this way the string network is 
constructed. 

Very recently there has been a trend to study the for- 
mation of defects and, in particular, strings, from a more 
field-theoretic viewpoint |||,^,||,||J^,^ and without adopt- 
ing a lattice Both efforts are very desirable and 
have shed further light on the problem. For example, the 
study of the evolution of fields during first order phase 
transitions (in two and three bubble collisions) has shown 
that the degree of freedom associated with can be 
important in the formation of defects. Generally, oscil- 
lations of |(/!)| lead to an excess of strings over a simple 
counting of the winding of a as in the VV algorithm. 
Also, it is often found that the number of defects is in 
excess of what the geodesic rule would imply. So, it seems 
that both assumptions of the VV algorithm might need 
to be relaxed. 

Another situation where the geodesic rule has been 
called into question is in the formation of gauge (local) 
strings In this case, the model in eq. (|]) is gauged. 
A consequence of this is that Va is no longer a gauge 
invariant quantity. For example, the contribution to the 
energy density no longer comes from the ordinary gradi- 
ent of a but from the covariant gradient: 

(Va - eAf 

where A is the gauge field. And this questions the 
geodesic rule that is based on the ordinary gradient of 
a. 

A significant advance has been made in Ref. where 
the collision of bubbles of true vacuum is studied analyti- 
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cally. The end result justifies the geodesic rule under cer- 
tain circumstances. However, in the analysis, the gauge 
field is not provided with any thermal properties of its 
own and is a passive player during the phase transition. 
In general, we would expect that the gauge field would 
have its own dynamics consistent with the thermal na- 
ture of the system. This could lead to the production of 
string in excess of what the geodesic rule would imply. 

FIG. 1. The filled circles show the string density as a func- 
tion of P in the Monte Carlo simulation. The curve shows the 

results of a direct calculation of the string density wficrc wc 
evafuated the expectated number of strings passing through a 
plaquette by summing over aif phase distributions and vaiues 
of n with appropriate probabiiity factors. 
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to be a random variable with the (Gaussian) probability 
distribution: 

/.ri+0.5 

Pn= dm 2v/^e-''(*"+2^™)' 

where, /3 G (0, oo) is a parameter. 

The form of P„ has two features which we think are 
desirable. First, for large (3, the algorithm reduces to 
the VV algorithm with the geodesic rule. Secondly, the 
Gaussian form is reminiscent of the Boltzmann suppres- 
sion factor with [3 being related to the inverse tempera- 
ture. On dimensional grounds, the connection of /3 and 
the temperature T is: 

T 

where ^ is the correlation length at the phase transition. 

FIG. 2. The fraction of infinite string as a function of the 

total string density. The fa-error bars arc found by running 
the simulation for each value of /3 f 00 times and then finding 
the spread in the fraction of infinite strings. The spread in 
values of the total string density is much smaller and has not 
been shown. 
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In this paper, we remedy the VV algorithm by relaxing 
the geodesic rule. Instead of calculating the gradient in 
the phase a on a lattice link by choosing the shortest 
path on the vacuum manifold, we choose the path using 
a probability distribution. More explicitly, if the phase 
at lattice site A is a a and that at lattice site B is as, 
then the difference in phase between sites B and A is: 

Aa = as — (XA + 27rn = 5a + 27rn 

where, n is any integer. The geodesic rule was that n 
should be chosen so as to minimize |Aa|. Here we take n 
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By changing the value of the parameter /?, we can now 
study the effects of relaxing the geodesic rule. This is 
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equivalent to allowing for the production of extra strings 
in bubble collisions, say due to fluctuations of |0|, or, due 
to gauge field fluctuations. 

We have performed the numerical simulations in this 
way on a cubic lattice together with the triangular dis- 
cretization of phases described in Q. Now it is possible 
for more than one string to pass through a plaquette in 
the lattice. This means that we can have a very high 
density of strings. Also, in connecting the string net- 
work, at every step there may be several choices to make. 
We make this choice randomly with equal probability as- 
signed to every possible connection. 

The total string density in the Monte Carlo simulation 
has been plotted as a function of the parameter f3 in Fig. 
1. On this figure we have also plotted the string density 
calculated by directly evaluating the expected number of 
strings passing through each plaquette. The two calcula- 
tions are in agreement, giving us confidence in the Monte 
Carlo simulations. 

FIG. 3. The loop length distribution for (3 = 0.05. 
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/? - lead to an excess of infinite string fraction. 

The graph in Fig. 2 can be fit by several different func- 
tional forms and the data does not convincingly point to 
any fit. In principal it should be possible to collect more 
data for higher values of the string density (lower values 
of f3) and then fit to a power law in the limit pmf / p 1. 
We have not been able to find the corresponding expo- 
nent since data collection rapidly becomes very expensive 
on computer resources as /3 is reduced. 

We have also plotted the number density of string loops 
as a function of the length (l) of the loop. For all values 
of /3, the result is the usual scale invariant distribution 
where the number density falls off as In Fig. 3 we 

show this length distribution (for /3 = 0.05) which is well 
fitted by: 
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To conclude, our results show that the implementation 
of the geodesic rule in the VV algorithm provides a lower 
bound on the fraction of infinite strings in lattice simu- 
lations. It also shows that any fluctuations in the gauge 
fields or magnitude of cj) that produce more strings, are 
likely to increase the fraction of inflnite string. 
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In Fig. 2, we show the string fraction in inflnite string 
{Pinf I p) as a function of the total string density p where 
an infinite string is defined to be one that is longer than 
the square of the lattice size. It is clear from the graph 
that departures from the geodesic rule - smaller values of 
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